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Abstract 

Analytical and numerical calculations show that a putative tempo- 
ral variation of the speed of light c, with the meaning of space-time 
structure constant cst, assumed to be linear over timescales of about 
one century, would induce a secular precession of the longitude of the 
pericenter w of a test particle orbiting a spherically symmetric body. 
By comparing such a predicted effect to the corrections Am to the 
usual Newtonian/Einsteinian perihelion precessions of the inner plan- 
ets of the Solar System, recently estimated by E.V. Pitjeva by fitting 
about one century of modern astronomical observations with the stan- 
dard dynamical force models of the EPM epehemerides, we obtained 
c/c = (0.5 ± 2) x 10~ 7 yr _1 . Moreover, the possibility that c/c ^ 
over the last century is ruled out at 3 — 12a level by taking the ratios of 
the perihelia for different pairs of planets. Our results are independent 
of any measurement of the variations of other fundamental constants 
which may be explained by a variation of c itself (with the meaning 
of electromagnetic constant cem)- It will be important to repeat such 
tests if and when other teams of astronomers will estimate their own 
corrections to the standard Newtonian/Einsteinian planetary perihe- 
lion precessions. 
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1 Introduction 

In this paper we will deal with the problem of effectively putting on the test 
an hypothetical time- variation of the speed of light c in a purely phenomeno- 
logical, model-independent way with local, Solar-System-scale astronomical 
observations. 

Varying Speed of Light (VSL) theories were proposed in recent times to 
accommodate certain features of the hot Big-Bang cosmology. In this sense, 
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the first modern VSL theory was put forth by Moffat in Ref. [I]; for other 
pioneering works see, e.g., Ref. [2] and Ref. [3]. Since then, this subject was 
dealt with by many authors investigating different aspects of it; see, e.g., 
Ref. [1] for an extensive review. Broadly speaking, such theories can be 
subdivided in two categories: those encompassing space-time variations of 
c, motivated by cosmology, and those where c varies with the energy scale, 
related to phenomenological quantum gravity. Subtle issues concerning fun- 
damental aspects of VSL theories have been recently discussed in Ref. [5], 
Ref. [6J and Ref. [Zj. 

From the observational point of view, the measured percent change [8, 
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— = (-7.2 ±1.8) x 1(T 6 (1) 
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of the fine structure constant^ 



a = f, (2) 
he 
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where q e is the rationalized electron charge and h is the Planck's constant, 
respectively, from quasar observations at redshift z ~ 0.5 — 3.5 was very 
important for VSL theories; indeed, the natural question arises: if a is 
varying, is such a change due to q e , h or c? By attributing Aa/a to a 
temporal variation of c, it follows 

Ac _ Aa 

c a 

According to the distinction of the many facets of c proposed in Ref. [5] , the 
c present here would be cem, i-e. the electromagnetic constant. Since 

Aa _ Opast ~ Qtoday < Q ,^ 
Q <^today 

the value of a was lower in the past; thus, c would have been larger in the 
past and it would be decreasing. By assuming a linear time dependence 

Ac c . 

— «~(t-*o <0, 5 
c c 

from eq. ([T]) it can be obtained 

- w (-8 ± 2) x 1(T 16 yr- 1 (6) 



For a review on the issue of the variation of a and other fundamental constants, see, 
e.g., Ref. p]. 
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for 



t - t « 9 Gyr 



(7) 



which approximately corresponds to the temporal interval spanned by the 
data analyzed in Ref. [10] , i.e. from 23% to 87% of the age of the universe. 
A tighter bound could be obtained from the constrain in the variation of a 
over the last 1.8 Gyr 

< 3 x 10~ 17 yr" 1 (8) 

from an analysis of the Oklo mine data |12j . 

Can locao (in space and time) astronomical observations tell us some- 
thing about the hypothesis that c undergoes temporal variations? The an- 
swer is, in principle, positive because the motion of the major bodies of the 
Solar System is governed by the dynamical equations of motion of classical 
general relativity in which c, playing the role of the space-time structure 
constant cst [S], is explicitly present; as we will see, a (slowly) time- varying 
c induces dynamical effects that can be tested with the latest planetary ob- 
servations independently of a. Of course, it must be borne in mind that 
such tests can only constrain c/c over timescales of about one century, cor- 
responding to the temporal interval covered by the modern astronomical 
observations of the major bodies of the Solar System which are used to 
construct the present-day highly accurate ephemerides. 



Q 
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2 The dynamical effects of c/c on the orbital mo- 
tion of a test particle 

We will follow a phenomenological approximation, without working in any 
specific VSL theoretical framework. By inserting eq. © into the 1PN 
gravitoelectric acceleration |14j of order 0(c~ 2 ) 



LlPN 



GM 

c 2 r 3 



/ AGM 



v 2 I r + 4(r • v)v 



(9) 



which causes the well-known Mercury's perihelion precession of 43.98 arcsec 
cy _1 , one gets 

(t-to) 



AA 



1PN 



LlPN: 



(10) 



2 For a strategy to combine local and cosmological tests of varying fundamentals con- 
stants like a and G see Ref. 1131. 
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here and in the following c = cq = c(to). Note that, according to the 
distinction of the many facets of c proposed in Ref. [5], the quantity varying 
here is the space-time structure constant c^t, which is, in principle, not 
related to the electromagnetic constant cem- It can be shown that the first 
term of eq. is proportional to c| T /c E , while the other two terms are 
proportional to l/c E , where ce the Einstein space-time matter constant [5]; 
however, in order to get the correct Newtonian limit for gravity [5], we will 
assume cst = ce- Our approach, which has the merit of making direct and 
unambiguous contact with the observations giving definite answers^], might 
be criticized from a theoretical point of view as, perhaps, too naive; indeed, 
as pointed out by Jordan |15l 116] . in general, it is not consistent to allow 
a constant to vary in an equation that has been derived from a variational 
principle under the hypothesis that this quantity is constant; one needs to go 
back to the Lagrangian and derive new equations with the constant treated 
as a dynamical field. However, whatever the temporal evolution of c(i) 
may be, the approximation of eq. © is adequate for the practical purpose 
of testing it over relatively short timescales like the last century in which 
modern astronomical planetary observations were collected. Incidentally, let 
us note that, in this case, certain observational issues 0E] can be neglected, 
at least from a practical point of view. Indeed, from 

dr = ^dt, (11) 
c 

where c is the space-time structure constant cst> by assuming a linear time 
variation of it, the following shift in the measured proper time would occur 
for a static field 

At c At . . 

r c 2 

over At = 100 yr and by assuming for cst the same rate of change of cem 
obtained from the Oklo natural reactor data for a of eq. ([8]), it turns out 



At 



1.5 x 10" i5 , (13) 

which is basically undetectable given the present-day accuracy in realizing 
the SI second by the Bureau International des Poids et Mesures (BIPM), 
i.e. pi] 3 x 10~ 15 . 

From a dynamical point of view, A Aipn can certainly be considered as a 
small perturbation with respect to the Newtonian monopole over timescales 



3 It should be recalled that the observational basis of VSL phenomenology is, at present, 
quite meager. 
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of about 100 yr, as it will be a posteriori confirmed by the bound on \c/c\ 
that we will obtain with such a hypothesis; the same holds if one uses eq. © 
derived from Aa/a. Thus, eq. (|10p can be treated perturbatively with the 
standard Gauss [18] approach which is valid for any perturbing acceleration, 
whatever its physical origin may be. In order to evaluate the orbital effects of 
a generic small disturbing acceleration W, it is customarily projected onto 
an orthonormal frame K co-moving with the test particle. The mutually 
orthogonal unit vectors r , f, t> of K pick out the radial, transverse and 
normal directions, respectively; f and t are in-plane with f directed along 
the particle's radius vector, while v is out-of-plane, directed along the orbital 
angular momentum. The Gauss equations for the variations of the Keplerian 
orbital elements are 



da 
~dl 
de 
If 
dJ_ 
It 
dtt 
~dt 
duo 
~dt 
dM 
~dt~ 
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where a, e, /, u and M are the semi-major axis, the eccentricity, the incli- 
nation, the longitude of the ascending node, the argument of pericentre and 
the mean anomaly of the orbit of the test particle, respectively. The angle / 
is the true anomaly reckoning the instantaneous position of the test particle 
along its orbit with respect to the pericentre, u = uj + f is the argument of 
latitude, p = a{\ — e 2 )is the semi-latus rectum and n = \J GM/a? is the 
un-perturbed Keplerian mean motion related to the un-perturbed Keplerian 
orbital period by P\y = 2n/n. For the following calculations it is more con- 
venient to use the eccentric anomarjQ E in terms of which the un-perturbed 
Keplerian ellipse at epoch to can be written as 



r = a(l — ecosE), 



(20) 



It is defined by M — E — e sin E. 
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Table 1: First row: eccentricities e of the inner planets of the Solar System. Second 
row: numerically calculated values of F(e) according to eq. (|26[) . 



Mercury 


Venus 


Earth 


Mars 


e 0.20563069 


0.00677323 


0.01671022 


0.09341233 


F(e) 527.063 


5862.73 


2554.07 


714.504 



cos E — e , . 

cos/ = - -, 21 

1 — e cos E 

. Vl - e 2 s'mE 

sm/ = — =-, 22 

1 — e cos E 

To obtain the secular, i.e. averaged over one orbital revolution, effects of 
W, it has to be evaluated onto the unperturbed Keplerian ellipse with the 
aid of eq. (|2T)j)-eq. ([22"]) and inserted into the right-hand-side of eq. (fT4"]) -eq. 
(|19P ; then, an integration with respect to t over an orbital period has to be 
performed by using 

dt n - eC ° sE ^dE. (23) 



In the case of eq. (fit)]) , with 



E — e sin E ... 

t-t = , 24 

n 



the Gauss equation for uj yields 



«-(l)(fr)^ 



with 



71^2 2(2, + e 2 )cosE + e{-lb + 7cos2E)(-E + esmE) 

Fie) = / — 3 dE; 

e J (1-ecosE) 3 

(26) 

in Tabled] we quote the numerically computed values of F{e) for the inner 
planets of the Solar System. Note that eq. (|25p holds also for the longitude 
of pericentre w = + uj; indeed, since W u = 0, from eq. (fTB]) and eq. (ITT]) 
turns out that ^/y = (tl^ = 0. Moreover, it is not possible to attribute 
the pericentre precession to a re-scaled time-varying gravitational constant 
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Figure 1: Numerically integrated trajectory of a fictitious planet around the Sun af- 
fected by the perturbing acceleration of eq. (|10|) in addition to the Newtonian monopole 
—GM/r 2 . A positive value large enough (c/c = 10 4 yr -1 ) to sufficiently enhance 
the pericentre precession has been chosen for c/c. The initial conditions chosen are 
xo = f min = a(l - e), j/o = 0, z = 0, x = 0, y = »m» = na A /(l + e)/(l - e), 
io = with a = 1 AU, e = 0.85: the motion of the planet along the orbit is anticlockwise. 
The temporal interval spanned by the integration is 10 yr. The retrograde (i.e. clockwise) 
precession, as predicted by eq. (|25p . is clearly visible. 

because, in this case, also the Newtonian monopole —GM/r 2 would be fic- 
titiously affected. Our analytical result for the pericentre precession of eq. 
(|25p is also qualitatively confirmed by a numerical integration of the equa- 
tions of motion of a fictitious planet around the Sun in which the magnitude 
of c/c has been purposely set to value large enough to visually inspect the 
rotation of the orbit in its plane, as depicted in Figure Q] 

Let us stress that eq. ([25]) is different from the precession obtained 
by Magueijo by investigating in Ref. [TU] spherically symmetric solutions 
to a definite covariant and Lorentz-invariant VSL theory [20]; indeed, the 
Magueijo's effect is equal to the usual 1PN precession^ multiplied by an 
adimensional factor, i.e. 

3nGM /4 6 2 \ 3(GM) 3 / 2 /4 6 2 \ 

WVSL " ~c 2 a(l-e 2 ) U « / "cW2(i_ e 2) 1,3 7 J ' (2?) 

where b and k are, in turn, numbers [20U19] . presumably of some fundamen- 
tal nature, accounting for the dynamical evolution of c. Simple dimensional 
considerations show, in fact, that eq. (|27p does not look like the formula 
one would reasonably expect for a weak-field dynamical precessional effect 
induced by a (slow) time variation of c. Indeed, the basic ingredients that 

5 The precession per orbit is shown in Ref. [19]; in order to compare it with our results, 
it must simply be divided by the Keplerian orbital period P b = 2n/n. 
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should intuitively enter such a formula are the two lengths GM/c 2 and a, 
characterizing the problem at hand, a possible adimensional function of the 
eccentricity e and a quantity Q having the dimensions of the reciprocal of 
time; [Q] = T _1 . Now, possible candidates for Q are the orbital frequency 
n and, of course, c/c which is the cause of the effect looked for; excluding 
quadratic terms in c/c, the most natural choice seems to be just Q = c/c. 
Stated differently, it would be possible to express eq. (f25|) as the standard 
1PN precession times an adimensional factor £, but the latter one would be 



where there is only one dimensional parameter related to the variation of 
c, i.e. its percent first derivative, while the other dimensional quantity, 
specific to the system considered, is the planet's orbital frequency. As we 
will see later, the dependence on a and e is crucial for the confrontation with 
observation-related quantities. 

3 The confrontation with the observations in the 
Solar System 

By suitably using the perihelia of the inner planets of the Solar System it 
is possible to constrain c/c over timescales of the order of about 1 century 
and even rule out the hypothesis that it may have a non-zero value, at least 
in the last century. 

The astronomer E.V. Pitjeva has recently fitted almost one century of 
planetary data of various types with the dynamical force models of the EPM 
ephemerides estimating, in the least-square sense, several parameters; the 
modelled dynamical features include [21] all the most relevant Newtonian 
effects (N-body mutual perturbations among the major bodies of the So- 
lar System, Sun's oblateness, 301 large asteroids, massive ring lying in the 
ecliptic plane accounting for the small asteroids) and the general relativistic 
Schwarzschild-like accelerations in the harmonic gauge. Among the various 
solutions obtained, in one of them she also phenomenologically estimated 
corrections \21\ [22] Azu to the standard Newtonian-Einsteinian precessions 
of the longitudes of the perihelia^) of the inner planets by keeping the usual 
PPN parameters fixed to their general relativistic values. By construction, 

6 Strictly speaking, the perihelia are not observables; they can be computed from the 
measured quantities which are ranges, range-rates and angles like right ascension and 
declination. 




(28) 
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Table 2: First row: estimated corrections Alb to the Newton/Einstein perihelion preces- 
sions of the inner planets, in 10 -4 arcsec cy -1 according to Table 3 of Ref. [21] (Mercury, 
Earth, Mars). The result for Venus has been obtained by recently processing radiometric 
data from Magellan spacecraft [52]. In square brackets we quote the formal, statistical 
errors resulting from the least-square estimation process. In the text we used the re- 
scaled errors. Second row: nominal general relativistic Lense-Thirring precessions zuur, 
10 -4 arcsec cy -1 . Such effects, not included in the dynamical force models of the EPM 
ephemerides and, thus, present, in principle, in the deterimed Aro, must be subtracted 
from the corrections Avj in order to single out the anomalous exotic effects induced by 
neither classical mechanics nor standard general relativity. 





Mercury 


Venus 


Earth 


Mars 


Am (10- 4 " cy- 1 ) 


-36 ± 50[42] 


-4±5[1] 


-2±4[1] 


1±5[1] 


w LT (10- 4 " cy" 1 ) 


-20 


-3 


-1 


-0.3 



such corrections Aw, shown in Tabled account, in principle, for any stan- 
dard, i.e. Newtonian and/or general relativistic, or exotic un- modelled/mis- 
modelled forces. In order to have the fully non-relativistic, exotic effects, 
the Lense-Thirring precessions, not modelled in the EPM ephemerides, have 
to be subtracted from the estimated corrections, i.e. one has to use 

Am* = Am - m LT . (29) 

Now, Am* can fruitfully be compared to eq. (125j) to constrain c/c by as- 
suming that they are entirely due to the putative dynamical effects due to 
the first derivative of c. 

By letting c/c be a free parameter, we can constrain it by comparing eq. 
(|25p and Tabled] to the estimated Am* quoted in Table [2] From a weighted 
mean of the values of c/c obtained with the four inner planets it turns out 

- = (0.5 ±2) x 10" 7 yr _1 , (30) 

compatible with eq. (jBJ). Our result, obtained without considering Aa/a 
and valid for the last century, is very conservative and pessimistic: indeed, 
we did not use the mere formal, statistical errors in Am and we linearly 
added the errors 5Am and Sm^/ C in constructing the total uncertainty in 
\Am* - m d / c \- 

By suitably combining the perihelia of various pairs of planets it is possi- 
ble to perform a more stringent test of the hypothesis that currently c/c ^ 0, 
independently of its origin and magnitude. Indeed, eq. ([25]) and eq. ([26]) 
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Table 3: First column: pair of planets A and B. Second column: observationally 
determined ratios II = Azu\/ Avj^, for A and B. Third column: theoretically predicted 
ratios A = Faob/-Fb«a f° r A and B. Fourth column: T = |II — ,4|/<5II; F > 1 means that 
II 7^ A within the errors. It turns out that the uncertainties in e and [23] a, are completely 
negligible in evaluating the errors in II — A which are, instead, dominated by 511. 



A B n AT 

Venus Mars -0.8 ± 6.8 17.3 3 

Earth Mercury 0.06 ± 0.44 1.87 4 

Venus Mercury 0.06 ± 0.50 5.95 12 



yield a function of a and e which represents a distinctive signature of the 
dynamical effects of c/c, irrespectively of its size; moreover, it is important 
to note that c/c enters eq. (j25l) as a multiplicative factor. Thus, by taking 
the ratios A of eq. (|25p for different pairs of planets A and B it is possible to 
construct theoretical predictions which are, at the same time, independent 
of the magnitude of c/c and still retain a pattern characteristic of c/c itself. 
Thus, A can be compared to n = Aw^/Aw^ for the same pairs of planets: 
if A 7^ II within the errors, i.e. if \A — H| ^ within the errors for some 
of the pairs considered, we must reject the possibility that c is nowadays 
varying according to c/c / 0. The results are in Table the hypothesis 
c/c / during about the last century must be rejected at more than 3 — a 
level. Also in this case, our test is conservative because we evaluated the 
uncertainty in II as 

It must be noted that, since Aw are observation-related quantities, it is 
perfectly meaningful to take their ratios II; the fact that 5 Am /Aw > 1 
simply means that Aw can still have a non-zero value smaller than 5 Aw 
and that II is compatible with zero within the errors. 

It maybe interesting to note that the perihelion precession of eq. (|27|) 
by Magueijo [19J would survive the test of the ratios of the perihelia. 

4 Discussion and conclusions 

In this paper we phenomenologically put on the test the hypothesis that 
the speed of light c, with the meaning of space-time structure constant 
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cst, can vary over timescales of about one century. We analytically worked 
out the dynamical effects induced by a linear variation in time of c on the 
motion of a test particle orbiting a spherically symmetric body finding that 
the longitude of pericentre w undergoes secular precessions; a numerical 
integration of the equations of motion qualitatively confirmed this result. As 
expected from simple dimensional considerations, the expression obtained 
for w is proportional to the product of c/c by (GM/c 2 a)F(e), where F(e) is 
a specific adimensional function of the eccentricity e. We compared such a 
theoretical prediction to the recently estimated corrections to the standard 
Newtonian/Einsteinian perihelion precessions for the inner planets of the 
Solar System, obtained by analyzing the last century of data, finding c/c = 
(0.5 ± 2) x 10~ 7 yr _1 . Moreover, by taking the ratios of the computed 
anomalous perihelion precessions for different pairs of planets we were able 
to obtain a prediction independent of c/c itself and still retaining a pattern 
characteristic of it. The confrontation of such predicted ratios with the 
ratios of the observationally determined corrections to the usual perihelion 
precessions ruled out the hypothesis that c/c / in the last century at 
3 — 12cr level. Our result is independent of any measured variations of other 
fundamental constants which could be related to a variation of c itself (with 
a different meaning like, e.g., that of electromagnetic constant cem)- If and 
when other teams of astronomers will estimate their own corrections to the 
standard perihelion precessions it will be possible to fruitfully repeat this 
tests. 
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